SOME OBSERVATIONS ON THE SIMPLEX 

IGOR RIVIN 

Abstract. We discuss some properties of the set of simplices in 

1. An inspirational example 

Consider triangles in the Euclidean plane E^. The set of congruence 
classes of triangles is parametrized by their three sidelengths (say a, 
b, c), while the set of similarity classes of triangles is parametrized 
naturally by the triples of angles a, /3, 7. 

What can be said about the possible values of the triple a, b, c? It is 
clearly necessary that the triangle inequalities be satisfied: 

a > 0, 6 > 0, c> 0, 
a < b + c, b < a + c, c < a + b. 

Since, given a triple T = (a, b, c) satisfying the three triangle inequali- 
ties we can easily construct the triangle with the sidelengths prescribed 
by T, it follows that, as parametrized by the lengths of sides, the set of 
non-degenerate triangles in is a convex cone. A similar result holds 
for the angles: a triple a, /5, 7 of positive numbers gives the angles of a 
non-degenerate triangle if and only if 

(1) a + /5 + 7 = 7r. 

One can ask whether results of this type extend to simplices of higher 
dimension, and we shall attempt to give as complete an answer as 
possible. 



2. Some negative results 

First, we shall see that things appear to become less simple. Firstly 

Theorem 2.1. Let Iq^i, ■ ■ ■ ,ln-i,n be positive numbers. In order for 
these to be the edge-lengths of a non-degenerate simplex vq, . . . ,Vn in 
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E", SO that lij = d{vi,Vj) it is necessary, but not sufficient that for any 
Vi, Vj, Vk the lengths kj, k^k, lj,k satisfy the requisite triangle inequalities. 

Proof. The necessity is obvious. To prove the insufficiency, we con- 
struct an example as follows: Let = 3, and let 

h,2 = ^1,3 = ^2,3 = I5 

while 

^0,1 = h,2 = ^0,3 = 1/2 + ^■ 

The reader will verify that for any e > all possible triangle inequal- 
ities are satifsfied, whereas for a sufficiently small e no simplex with 
prescribed edgelengths exists. □ 

Things are even worse than that: 

Theorem 2.2. Let lij, < i, j < n be as in the statement of The- 
orem Wl[ The set of simplices in E" (for n > 2) is not convex when 
parametrized by the Uj. 

Proof. (The example is due to Peter Frankel of Budapest). Consider 
two simplices A and B as follows: 

Vi(^) = /o,2(A) = loAA) = loAB) = loAB) = loAB) = 1, 

while 





= e, 




= V2, 




= V2, 




= V2, 


12AB) 


= V2, 


hAB) 


= e. 



It is a simple geometric exercise to show that A and B really and truly 
exist (for any e > 0) while, for e small enough, there is no simplex 
C with lij{C) = kj{A) + li,j{B), VO < i < j < 3. (alternatively, 
the reader can peek ahead for a hint to a non-geometric proof of non- 
existence of C.) □ 

We shall have to resort to less geometric methods to see what is 
really going on. 
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3. Linear algebra to the rescue 

Let A be a simplex with vertices = 0, vi, . . . , in E". The volume 
of the simplex A is non-zero if and only if the vectors Vi, . . . ,Vn are 
linearly independent; in fact, its volume is given by 

(2) volA = ^ dct V, 

nl 

where V is the matrix whose columns are the vectors Vi, . . . ,Vn. 
We now write the so-called Gram matrix of A: 

G{A) = V'V. 

Equivalently, 

Gij{A) = {vi,Vj). 

Theorem 3.1. An nxn matrix G is the Gram matrix of a simplex A 
in E" if and only if G is symmetric and positive definite. 

Proof. First we prove the "if" direction. The symmetry of the inner 
product in E" imphes that G is symmetric (alternately: (V^Vf — V^V.) 
And for any vector x e E": 

^^^(A)^ = x^V^Vx = {Vx, Vx) = \\Vx\\ > 0, 

if V is non-singular and x is non-zero, hence positive definiteness. 

To show the "only if" direction, let G be a positive-definite sym- 
metric metrix. That implies that there is an orthogonal matrix S and 
a diagonal matrix D all of whose diagonal elements are positive, such 
that G = SWS. Write F = VD (the elements of F will be simply the 
non- negative square roots of the corresponding elements of D) . Then, 
it is easy to see that G = S'DS - S^F^FS = {FSy{FS). Setting 
V = FS wc obtain our simplex. □ 

The fact, used in the proof of the above theorem, that a symmetric 
matrix has an orthogonal matrix of eigenvectors is essentially equivalent 
to the following 

Theorem 3.2 (Rayleigh-Ritz characterization). The smallest eigen- 
value Ao of a symmetric matrix A can be written as 

An = min (Ax, x). 
11^11=1 

from which we have 
Corollary 3.3. 

Xo{A + B)>Xo{A) + Xo{B) 

( in other words, Aq is a concave function on the set of symmetric ma- 
trices. 
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Proof. By the Rayleigh-Ritz characterization, 

Xo{A + B) = mill {{A + B)x,x). 

\\x\\=l 

Denote the unit vector achieving the minimum by xa^b- Since 
{Axa+b, xa+b) > min {Ax, x), 

\\x\\=l 

and similarly 

{Bxa+b,xa+b) > min (5a;, x), 

|x||=l 

the assertion follows. □ 

And the final corollary: 

Corollary 3.4. The set of positive- definite symmetric matrices is a 
convex cone in the space ofnxn symmetric matrices (which can be be 
naturally identified with Mv 2 J. 

Proof. This is an immediate consequence of Corollarv 13.31 since a ma- 
trix is positive-definite if and only if all of its eigenvalues are posi- 
tive. □ 

Corollarv 13. 41 can be restated as saying that the set of Gram matrices 
of simplices is a convex cone. This does not see very useful at the 
moment, since the entries of the Gram matrices are some obscure scalar 
products, but this can be easily rectified by observing the 

Polarization identity: 

{v,w) = - {\\v\\'^ + II w IP — ||f — u'lp) . 
When applied to the Gram matrix, this shows that 

where are the edge lengths as before, and U^i = for any i. We 
see that the Gram matrix is a linear (matrix valued) function of the 
squares of the edge lengths of the simplex A, and hence we can put all 
of the above together to see that: 

Theorem 3.5. The set of non- degenerate simplices in E" is a convex 
cone when parametrized by the squares of the edgelengths. 

We have now recovered our convexity of the set of simplices, which, 
in view of the gloomy Section 2 can already be viewed as a success. 
But, as they say. Wait! There is MORE! 
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Remark 3.6. The following is related to the results of this section, and 
follows from a beautiful theorem of I. J. Schoenberg ( |Schoenberg37| ) : 
the set of simplices parametrized by the logs of the edge lengths is 
starshaped with respect to the origin (corresponding to the regular 
simplex) . 

4. Volume 

In this section we shall use the following simple (as the reader will 
see momentarily) result: 

Theorem 4.1. The function logdetA is a concave function on the 
cone of positive definite symmetric matrices. 

Proof. Concavity is equivalent to concavity on all lines, which is in turn 
equivalent to the statement that 

d^\ogdet{A + tB) 

(3) 



dt^ 



<0, 

i=0 



for any positive definite A and any symmetric B. To show the inequality 
Q we have to be able to differentiate the determinant (or its logarithm 
directly). A couple of ways of doing that will be shown below, but for 
now we will accept the following as a fact of life: 

(4) ^^^^^^^^±^ = triiA-,tBrB). 

From which it follows easily (by linearity of trace) that 
\ogdet{A + tB) 



(5) 



dt^ 



= -tT{A-^BA-^B). 

t=o 

Now, since determinant is invariant under conjugation, we can assume 
that A is diagonal: A = diag(Ai, . . . , A„). Under that assumption, di- 
rect computation shows that 

A-'B]= V 



tTiA-'BA-'B)= J2 



0<i,j<n * 

Since all the summands are nonnegative, and at least one is positive 
(since B 0), the theorem is proved. □ 

We finally obtain the following: 

Theorem 4.2. The volume of any fixed given face (of any dimen- 
sion) of simplices parametrized by the squares of the edgelengths is log- 
concave. 
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Proof. The result for the top-dimensional face follows from the expres- 
sion (0) and Theorem 14. II above. Since the Gram matrix of any given 
face is a submatrix of the Gram matrix of the whole simplex, the result 
follows. □ 

Using the usual inequality theory (see fG ardner02] for details) we 
obtain the following result of a Brunn-Minkowsky type: 

Corollary 4.3. The n-th root of the volume of n- dimensional simplices 
parametrized by the squares of the edgelengths is a concave function. 

These results can be used to produce hitherto unknown characteri- 
zations of the regular simplex. First, an observation: 

Lemma 4.4. Let Q be a compact convex set in an affine space, G a 
group of automorphisms acting on Q, and f a concave function invari- 
ant under G (that is, f{G{x)) = f{x), \/x e fl). Let y be a point 
where f achieves its maximum. Then y is invariant under G ( that is, 
g{y) = y, \fgeG.) 

Proof. Suppose that there is a g E G such that z = g{y) ^ y. Since VL 
is convex, the segment S = [y, z] is contained in Q. The function / is 
concave on S, and so 

/ (^) > ^ f(y), 

contradicting the assumption that y is a maximum of /. □ 
And now, the applications: 

Theorem 4.5. Let Aa be the set of n-dimensional simplices, such that 
the sum of the squares of their edgelengths is equal to a. Let 

Pk= n ^(/)- 

k- dimensional faces f of A E Aa 

Then Pk is maximized when A is a regular simplex (with edgelengths 
2a/{n + l){n). 

Theorem 4.6. Let Aa be the set of n-dimensional simplices, such that 
the sum of the squares of their edgelengths is equal to a. Let 

Sk= Yl ^'^'(/)- 

k- dimensional faces f of A E Aa 

Then Sk is maximized when A is a regular simplex (with edgelengths 
2a/{n + l){n). 

Proof. The Theorems above follow immediately from Theorem l4 . 21 ( Gorol- 
larv 14. 3t respectively) and Lemma f4. 41 □ 
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5. How TO DIFFERENTIATE THE DETERMINANT 

In this section we demonstrate the truth of the formula (jH). Consider 
the family 

M{t) = A + tB. 

Clearly, 

det M{t) = det A det [I + tA-^B). 

The logarithmic derivative of det M{t) at is thus equal to the deriv- 
ative of det(/ + tA~^B) at 0. To simplify notation, write C = A~^B. 
The i-th column of I + tC equals ej + tCj, where Cj is the ith standard 
basis vector, and Cj is the ith column of C. Using the multilinearity of 
the determinant, we see that 

det{I + tC) = 
det(ei +tci,...,e„ + tc„) = 

n 

l + tJ2 det(ei, . . . , Ci, . . . , e„) + 0{t^) = 

i=l 

l + tr(C) + 0(t2), 

showing that 

ddet{l + tC) _ 
dt - 

which completes the proof. 

6. Random facts about matrices 

Definition 6.1. Let M be a matrix. The adjugate M of M is the matrix 
of cofactors of M. That is, Mij = (-1)'+^ det M'^ , where M'^ is M with 
the i-th row and j-th column removed. 

The reason for this definition is 

Theorem 6.2 (Cramer's rule). For any n x n matrix M (over any 
commutative ring) 

MM = MM = (det M)I{n), 
where I{n) is the n x n identity matrix. 
We also need 

Definition 6.3. The outer product of vectors v = [vi, . . . , f„) and w = 
{wi, . . . , Wn) is the matrix v ® w, defined as follows: 

(f w)ij = ViWj. 
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Remark 6.4. As the notation suggests, the outer product is actually a 
tensor prcduct, though it would be more correct to write v (8) w*. The 
Dirac notation for the outer product would be 1^)11(^1 while the Dirac 
notation for the inner product would be {v\w), this possibly explaining 
the name outer product. 

Consider an arbitrary vector x — {xi, . . . , We see that 

n 

(6) [{v w)x]^ = y^^VkWjXi = Vk{w,x), 

i=l 

so that 

(7) {v ® w)x = {w,x)v. 

We see that v w is a multiple of the projection operator onto the 
subspace spanned by v. In particular, in the case when \\v\\ = 1, the 
operator v<^v is the orthogonal projection operator onto the subspace 
spanned by v. Since v ® w is a rank 1 operators all of its eigenvalues 
are equal to 0. The one (potentially) nonzero eigenvalue equals {v,w). 
We now show: 

Theorem 6.5. Suppose that M has nullity 1, and the null space of M 
is spanned by the vector v, while the null space of is spanned by the 
vector w. Then 

M — cw <S> V, 

Proof. Since M is singular, we know that detM = 0, and so every 
column of M is in the null-space of M. so, letting rrij denote the ith 
column of M, we see that 

m, = diVi. 

However, M* = (^M^ so performing the computation on transposes 
we see that 

m*j = CiWi. 

We see that 

Mij = diVj = CjWi. 

Writing di = giWi, and Cj = hjVj, we see that, for every pair i,j, 
QiWiVj = hjWiVj. The conclusion follows. □ 

Theorem 6.6. The constant c in the statement of the last theorem 
equals the product of the nonzero eigenvalues of M divided by the inner 
product of V and w. 
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Proof. By considering the characteristic polynomial of M we see that 
the product of the nonzero eigenvalues of M equals the sum of the 
principal n — 1 minors. On the other hand, the principal minors of M 
equal the diagonal elements of M, so 

n n— 1 

i=l i=l 

□ 

Remark 6.7. By the discussion following Eq. [7| the product of nonzero 
eigenvalues of M equals 

det {M + w^v) 
{v,w) 

7. Back to simplices 

First we define the dual Gram matrix of a simplex S. First, let let 
fi be the unit outer normal to the i-th face. Then 

Gij = {fi, fj)- 

In other words, the ij-th entry of the dual Gram matrix is the cosine 
of the exterior dihedral angle between the i-th and the j-th face. 

Lemma 7.1. The dual Gram matrix G*{A) of a Euclidean simplex A 
is symmetric and positive semi- definite, with exactly one eigenvalue. 

Proof. The proof proceeds exactly as the proof of Theorem 13. II Any n 
of the vectors /o, . . . , /„ are linearly independent, and the correspond- 
ing nxn submatrix of the Gram matrix G* is positive definite, whence 
the result. □ 

The above Lemma makes one wonder what the null space of G* 
might be. Happily, there is a complete answer, as follows: 

Theorem 7.2. The null space of G* is generated by the vector A = 
{Aq, . . . , An), where Ai is the area of the face fi. 

Proof. This follows from the divergence theorem, which states (in the 
polyhedral case) that 

n 

(8) 5^A/. = o. 

i=0 

Taking the dot product of eq. ^ with fj, we see that 

n 

(9) Y.MMJ)) = o, 

i=0 



10 



IGOR RIVIN 



which is to say that the j-th coordinate of G*A vanishes. □ 
Corollary 7.3. With notation as above, 

A] %^ 

Proof. Immediate from Theorem 16.51 and Theorem 17.21 □ 

Remark 7.4. This result has now been generalized to simplices in 
by the Novosibirsk State University Geometry seminar. 
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